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We investigate the linear growth rate of cosmological matter density perturbations in a viable f{R) 
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I. INTRODUCTION 

Cosmological observations of distant la supernovae discovered that our Universe is undergoing an accelerated 
expansion period pj, [2]] , which is supported by other observations of the cosmic microwave background anisotropics 
and the large scale structure of galaxies [^-Q. These observations are explained by the cosmological model with the 
cosmological constant A. The cosmological constant can be regarded as the vacuum energy, however, the smallness of 
the observed value raises a fine-tuning problem [6j . To explain the cosmic accelerated expansion, many dark energy 
models have been proposed (see e.g., [3, Q and the references therein). Modification of the gravity theory is an 
alternative approach, for example, HR) model and the Dvali, Gabadadze, and Porrati (DGP) model in the 

context of the braneworld scenario |13j . 

Many authors have studied dynamical dark energy models p^ - [l7j . Dynamical dark energy models may have similar 
expansion rates to models of modified gravity, because modification of the gravity theory may affect the background 
expansion history. Therefore, the observations of the background expansion history alone are unable to distinguish 
between modified gravity and dynamical dark energy. The key to distinguish between modified gravity and dynamical 
dark energy is the growth of cosmological perturbations |18l - [23j . The growth history of cosmological perturbations 
can be tested with the large scale structure in the Universe. Many projects of large survey of galaxies are in progress 
or planned [23 - [33j , and these surveys might give us a hint in exploring the origin of the accelerated expansion of the 
Universe and the nature of gravity [34-47]. 

Cosmological perturbations in modified gravity models have been investigated by many authors [48l - l62| . In the 
present paper, we investigate the growth history of matter density perturbations in /(i?) models, /(i?) model is a 
modified gravity model, constructed by replacing the gravitational Lagrangian with a general function of the Ricci 
scalar R. The viable f{R) models have been proposed [6314671 ] . which explain the late-time accelerate expansion of the 
background Universe, and satisfy the local gravity constraints. The viable model which also explains an inflationary 
epoch in the early Universe is extensivel y prop osed [68l - l70j . For the local gravity constraints, the chameleon mechanism 
is supposed to play an important role [7ll473l]. By this mechanism, a field that modifies the gravity is hidden in the 
local region with high density. We note that a problem of the theory in the strong gravity regime is under debate 
[tI . 111]. Though the evolution of cosmological perturbations in f{R) models has been studied so far (76l - [85j . our 
investigation is focused on a new description of the growth rate for the f{R) model. 

This paper is organized as follows. In Sec. II, we briefiy review the viable f{R) models. In Sec. HI, we investigate 
the evolution of density perturbations in the f{R) model both numerically and analytically. We find that the growth 
rate of density perturbations can be characterized by a simple analytic formula, which approximately describes the 
growth rate in the scalar-tensor regime. The growth rate in the general-relativity regime is also investigated. In Sec. I 
V, we investigate a future prospect of constraining the f[R) model assuming a future large survey of weak lensing 
statistics on the basis of the Fisher matrix analysis. Section V is devoted to summary and conclusions. Throughout 
the paper, we use the unit in which the speed of light equals 1 and h= 1. 

II. A BRIEF REVIEW OF f{R) MODEL 
We briefly review f{R} model, which is defined by the action. 
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where G is the gravitational constant, and L^™) is the matter Lagrangian density. We consider the viable models, 
proposed in Refs. (63l - [67j . The viable models have an asymptotic formula at the late-time Universe {R ^ i?c), which 
can be written as 



f{R) = -\R, 



1 



R£ ^" 
R 



(2) 



where R^. is a positive constant whose value is the same order as that of the present Ricci scalar, and A is a nondimen- 
sional constant. Because the term XR^ plays a role of the cosmological constant, we may write Ai?c = 6(1 — Vlf))HQ, 
where -ffo is the Hubble constant and fJo is the matter density parameter. Note that we assume the spatially flat 
Universe. 

It is well known that = df{R)/dR plays a roll of a new degree of freedom, which behaves like a scalar field with 
the mass 

2 + 



where we defined fjiji = d?f{R)/dR^. Assuming <C 1 and R/rr ^ 1 for the viable model, the mass is simply 
m^ = l/{3fRR). 

We focus on the evolution of matter density perturbations in the f{R) model, whose Fourier coefficients obey (e.g., 

M) 

S + 2HS - AnG^spS = 0, (4) 

where the dot denotes the differentiation with respect to the cosmic time, H — d/a is the Hubble parameter, p is the 
matter mean density, and Gcs is the effective gravitational constant, which is written as 

Geff ^ 1 P/a^ 
G 3P/a2 + l/(3/flfl)' 

where k is the wave number, and a is the scale factor normalized to unity at present epoch (cf. [s^l). As is noted 
in the above, the physical meaning of — l/(3/i{fl) is the square of the mass of the new degree of freedom which 
modifies the gravity force. We have the general-relativity regime, Gcfr ~ G, for k/a ^ m, and the scalar-tensor 
regime, Goff = 4G/3, for k/a ^ m, respectively. Thus, the evolution of matter density perturbations depends on the 
wavenumber k, whose behavior is determined by the mass rri^ — l/(3fRR). 

For the Einstein de Sitter universe, the exact solution of Eq. Q is found in the literature [8l|. However, we consider 
the low density universe, where the solution of Eq. ((4]) is described in a different form in comparison with that of Isij . 
From Eq. ([2|), we have 

d^ f ( n\ n^^^ 



Furthermore, using the formulas XRc = 6(1 - no)H^ and R = [^o/a^ + 4(1 - f^o)] , we have 

(7) 



1 _ noH^ f>^Vf 1 4(l-r!o)^'"+' 



3/i?j?, 4n{2n + l)\2j \1 - J flo 
Denoting the wavenumber corresponding to the Compton wavelength 1 /m at the present epoch by kc , 

npHj (xV-( Qo y"+7, , 4(l-f^o) y"+^ 
^'^-4n(2n + l) [2 J [T^o ) I + J ' 

Equation ([7]) is rewritten as 

J__,2 / »oa-^ + 4(l-»o) \ 

sfRR no + Ail -no) J ■ 

We denote the growth factor by Di{a, k), which is the solution of Eq. Q normalized so as to be Di{a, fc) ~ a at 
a <C 1. The growth rate is defined by 

dlogDi(a,fc) 

/(a, k) = . (10) 

d log a 



3 



Using the growth rate /(a, k), Eq. Q is rephrased as 

where fljn{a) is defined by i7„j(a) — H^Qoa^^ /H^. We assume that the background expansion is weU approximated 
by the ACDM model, where the Hubble parameter satisfies 

A = -^nM (12) 

and the energy conservation equation 

din a 

Using Eqs. ([n|)-(IT3l) yields 



-3n^{a)il-nm{,a)). (13) 



which is useful to find an approximate solution, as we see in the next section. 



III. GROWTH OF DENSITY PERTURBATIONS IN f{R) MODEL 

In this section, we investigate the evolution of matter density perturbations in the f{R) model. In Sec. M A, we 
consider the scalar-tensor regime, /c/a 3> m, in which the wavelength is shorter than the Compton wavelength. In 
Sec. in B, we consider the general-relativity regime, k/a <^ m, in which the wavelength is larger than the Compton 
wavelength. 



A. scalar-tensor regime 



In the scalar-tensor regime, k/a ^ m, the effective gravitational constant becomes Gcff = 4G/3. In this case, we 
find that Eq. has the solution expressed in the form j86j 



/(a,fc) = /ol)™(a)^('^). 



where /o obeys /q + /o/2 = 2, therefore /o = (-1 + v33)/4, and 

7(a) =ECK1-^™ («))'. 

where Q is the expansion coefficients. The first few terms of 7(a) are 

9-V33 93-17^33 



7(a) 



(l-r!„,(a)) + 0((l-a„(a))') 



(15) 



(16) 



(17) 



6 666 

This can be generalized to the case when Gcs/G{= C) is a- constant value, in which the solution of Eq. ([Tl]) has the 
same formula as that of p^)) but with /o = (— 1 + \/l + 24^)/4 and 



7(a) = 



-41 + 24g + yi + 2Ai 
-70 + 48^ 



1 



-143 + 240(-35 + 24^2 



X (-41 + 24^+ Vl + 24e)(347- 17Vr+24e + 24e(-13+ Vl + 24e)) 
X [l - n„,{a)) + O {{I - n„^{a)f) . 



(18) 



Here we assume that Gcs/G{— ^) is constant, but we utilize this formula by replacing ^ with the right-hand-side of 
Eq. ©. 
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FIG. 1: Left panel: f(a,k) as a function of z(= 1/a — 1) for the model n = 1. The solid curves are obtained by solving the 
differential equation numerically, for K{= k/kc) = 10^, 10, and 1, respectively, from the top to the bottom. The long 
dashed curve, the dotted curve and the short dashed curve adopt the approximate formula, for K{— k/kc) ~ lO'^, 10, and 1, 
respectively, from the top to the bottom. Right panel: The growth index 7(0, fc) as a function of z[= 1/a — 1), corresponding 
to the left panel. The parameter of the curves is the same as that of the left panel. The curves correspond to K{= k/kc) ~ 
and 10, 1 respectively, from the bottom to the top. 



The left panel of Fig. [T] shows the growth rate /(a, k) as a function of the redshift z{— 1/a — 1). Here we adopted 
n = 1. The solid curves are obtained by solving Eq. ([TT]) numerically, for K{= k/kc) — 10^, 10, 1, from the top 
to the bottom, respectively. Here we assumed the background expansion of the Universe is the ACDM model with 
ilo = 0.28. The dot-dashed curve, the dotted curve, and the short dashed curve are the approximate formula up 
to 1st order of (1 — 51.,„(a)), for the wavenumber k/kc) = 10^, 10, 1, respectively. In the computation of the 
approximate formula, we adopted the right-hand-side of Eq. ^ as ^. One can see that the approximate solution 
approaches the exact solution at the late time of the redshift. 

Following the previous works (e.g., see (5^), the growth index 7(0, k) is introduced by 

/(a,fc) = r!„(ar('^^'=), (19) 

which is related with 7(0, k) by 

^("'^) = r-^4T+'^(«'^)- (20) 

In il,n[a) 

The behavior of the growth index j(a,k) in the scalar-tensor regime is well approximated by Eq. (j20p . as is 
demonstrated in the right panel of Fig. [U which plots 7(0, fc) as a function of the redshift z, for the wavenum- 
ber K{= k/kc) — 10^, 10, 1, from the bottom to the top, respectively. The solid curves are obtained by solving 
Eqs. ([TTT) and numerically, for K{— k/kc) = 10^, 10, 1, from the bottom to the top, respectively. The dot- 
dashed curve, the dotted curve, and the short dashed curve are the approximate solution of 7(0, fc) and ((20)) . for 
K{= k/kc) = 10^, 10, 1, respectively. One can see that the approximate solution approaches the exact solution at 
the late time of the redshift, however, the validity is limited to the late time of the small redshift. 

Let us discuss the valid region of the approximate solution. The left panel of Fig. ^ plots the redshift as a 
function of K{= k/kc) for n = 1/2, 1, 2, 3, 4, and 5, respectively, from the top to the bottom, where z^ is defined by 
the redshift when the difference of the growth rate becomes /('^pp'') — j(exac) _ q q3_ jjere /l^^'^) is the exact solution 
obtained by solving Eq. (ITU) numerically, while /(^pp'') is the approximate solution. Thus, the approximate solution 
of the growth rate approaches the exact solution after the redshift z^, which depends on k/kc as well as n. As n is 
larger or k/kc is smaller, Zx becomes smaller. For the case n > 4 and smaller value of k/kc, we have no solution of 
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FIG. 2: Left panel: Redshift when the difference of the growth rate becomes /(^'pp'') - = 0.03, as a function of 

K{— k/kc). The curves are n = 1/2, 1, 2, 3, 4, and 5, respectively, from the top to the bottom. Right panel: Transition 
redshift Zc as a function of K{= k/kc)- From the top to the bottom, curves are n — 1/2, 1, 2, 3, 4, and 5, respectively. Here 
we adopted the background expansion of the Universe is the ACDM model with Slo ~ 0.28. 



The above behavior is related with the transition redshift Zc, when the scalar-tensor regime starts, which we defined 
by fc(l + Zc) — m, i.e., 

The right panel of Fig.[2]plots Zc as function of K{— k/kc) for n — 1/2, 1, 2, 3, 4, 5, respectively, from the top to the 
bottom. Figure [2] shows Zx < Zc- Thus the approximate formula approaches the exact solution after the scalar-tensor 
regime starts. For the model with larger value of n, the Compton scale evolves rapidly. Then, the transition redshift 
Zc becomes small as n becomes large. For the smaller value of K(= k/kc), the transition redshift Zc becomes smaller. 
This is the reason why z^ is smaller, as n is larger or k/kc is smaller. Therefore, for the case when n is large and 
k/kc is smaller, the redshift when the approximate formula starts to work becomes later. For the case n ^ 2, the 
late-time behavior of the growth rate can be approximated by the approximate formula as long as K ^ 1. 

We here mention the relation between the parameter kc and the parameter fjio adopted in Refs. [g^, H^l, in 
which the case n = 1/2 is investigated. In this case, |/i?o| = 2(1 — 3ilo/4)HQ /k^. For |/ho| — 10~^ — 10~^, we have 
kc — 0.04 — 0.4 hMpc^^. The scalar-tensor regime appears rather earlier in this model, as shown in Fig. 



B. general-relativity regime 



In this subsection, we consider the growth rate of density perturbations at the early time epoch of the Universe, 
a <C 1, adopting the approximation. 



3fRR An{2n+l)\2j \l - Qq J \a 
which yields the simple form of the effective gravitational constant 



(22) 
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where 



N = 2n 



4 

+ 3' 



An{2n + 1) 



2n 



no 



2n+l 



(24) 
(25) 



As mentioned in the previous section, 3/i?,ij has the meaning of the square of the Compton wavelength. Thus this 
model can be regarded as the model that the Compton wavelength simply evolves as 1/m = a'^^/^+^/fco- 

In the case when is a positive integer, we derive an approximate solution of Eq. in an analytic manner. 
With the use of |T9|), one can rewrite Eq. (fT4|) as 



3Sl™(a)(l - ri„(a))lnfi„(a)- 



1 



3 Go 



1-7 _ 



37 + 2 = 0. 



In a straightforward manner, we find the solution for 7(0, k) expanded in terms of 1 — r2m(fl)j as follows. 

7(a,A:) = ^CHfl,fc)(l-^^m(a))^ 



(26) 



(27) 



where (a, k) is the expansion coefficient. For example, for TV = 1, we find 



7(a, fc) 



6 K^no 
IT " 11(1 - fio) 

15 ISlK^flo 



2057 

where K = k/k^. We also have 



4114(1 -f7o) 4114(1 -f7o)2 



[l ~ n,n{a)) + O {{I - n^{a)f 



(28) 



15 



2057 17(1 - n^f 
4205 643X2^72 



1040842 8602(1 - r^o) 



(1 - n„,{a)) 

(l-l]™(a))2+0((l-a„,(a))3). 



(29) 



for N ■ 



6 

IT 



(l-a„(a)) 



15 
2057 
4205 
1040842 ~ 23(1 



(30) 



for iV = 3, 



l{a,k) = -7 + 



15 



2057 
31449595 



(l-l]„(a)) + 



ii-n^{a)f 



11288972332 29(1 - f^o)'' 



4205 
1040842 ^ 

(l-f}„(a))3 + 0((l-f]„(a))4), 



(31) 



for = 4, respectively. 

Figure [3] demonstrates the validity of the approximate formulas, by plotting the relative difference between the exact 
solution / and the approximate solution fij^ as a function of the redshift. The four panels assume A^ = 1, 2, 3, 4, 
respectively. In each panel, the cases of the wavenumber K{— k/kc) = 1, 0.1, 0.01, are plotted. The solid curve, the 
dotted curve and the dashed curve correspond to AT = 1, 0.1, 0.01, respectively. For A^ = 1, we used the approximate 
formula up to the 1st order of (1 — rim(a)). For A = 2, we used the approximate formula (|29l) up to the 2nd 
order of (1 — il„i(a)). For A^ = 3, we used the approximate formula (pO| up to the 2nd order of (1 — ri„i(a)). For 
A = 4, we used the approximate formula ([31]) up to the 3rd order of (1 — rim(a)). Even if we adopted higher order 
term of (|28p ~ pip . the approximate formula only slightly improves the accuracy for AT = 1. The approximate formula 
is valid for K ^ 0.1. 
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FIG. 3: The relative difference of the growth factor between the exact solution /, which is obtained by solving Eq. 
numerically, and the approximate solution in the form QZi, as a function of z. The four panels correspond to TV = 1, 2, 3, and 
4, respectively. In each panel, the curves correspond to A" = 0.01, 0.1, 1, respectively, from the bottom to the top. 



IV. CONSTRAINT ON f{R) MODEL FROM WEAK LENSING SURVEY 



Cosmological constraints on the f{R) model have been investigated in Refs. [6J,|87h89|. The weak lensing statistics 
is useful to obtain a constraint on the growth history of cosmological density perturbations observationally. We now 
consider a prospect of constraining the f{R) model with a future large survey of the weak lensing. To this end, we 
adopt the Fisher matrix analysis, which is frequently used for estimating minimal attainable constraint on the model 
parameters. To be self-contained, we summarized the fisher matrix analysis in the Appendix (see also [lO^, and the 
references therein). Here we focus on the constraint on the Compton wavenumber parameter kc defined by Eq. (|8]) 
or (|9l). In this analysis, we obtained the growth rate and the growth factor by numerically solving Eq. and 



Di(a, k) — aexp 



r^(/(a',fc)-l) 
Jo a 



(32) 



without using the approximate formula. 

We briefly review how the signal of the weak lensing reflects the modification of the gravity in the f{R) model. In 
the Newtonian gauge, the metric perturbations of the Universe can be describe by the curvature perturbation $ and 
the potential perturbation , 



ds^ = - (1 + 2*) dt^ + a^(t) (1 + 2$) dx^ 



(33) 



In the f{R) model, the relations between the two metric potentials and the matter density perturbations are altered. 
In the subhorizon limit, the f{R) model yields (e.g., [90] and references therein) 



fc^v]/ = -4:TrG^{a,k)a'^pS 
- = -i^{a,k), 



(34) 
(35) 
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FIG. 4: Left panel: The 1-sigma contour in the (n — A) plane. The linear modeling for Pmass{k, z) in the range of 10 < 
/ < 10"^ is used. The target modes are (n, A) = (3, 4), (2, 2), (5, 2), and (2,8), respectively. The other target parameters are 
Wo ~ —1, Wa = 0, flo ~ 0.28, Qb ~ 0.044, h — 0.7, as = 0.8, and Us — 0.95. The solid curve corresponds to kc = 0.2feMpc~^. 
Right panel: Same as the left panel, but with the nonlinear modeling for Pmass(fc, z) of the range of 10 < Z < 3 x 10''. 



with 

na,k) = „ — rj, 37) 

where we used \fji\ <^ 1 and RJrb, ^ 1. Equation ([M)) is the modified Poisson equation. In general relativity, 
^ = u=l. With Eqs. ([Ml) and jUj, we have 

fc2($_^') = 87rGa2p,5. (38) 

Thus, this relation between $ — and S is the same as that of the general relativity. The signal of the weak lensing 
is determined by $ — ^' along the path of a light ray. Therefore, we only consider the effect of the modified gravity 
on the matter density perturbations of Eq. (|4]) for elaborating the weak lensing statistics. 

In the present paper, the modified gravity of the f{R) model is supposed to be characterized by n and fcc(or A). 
We perform the Fisher matrix analysis with the 9 parameters, n, A(or fcc), wq, Wa, fib, h, A, and n^, where J7f, 
is the baryon density parameter, Us is the initial spectral index, A is the amplitude of power spectrum, wq and Wa 
characterize the background expansion history and the distance-redshift relation [see Eq. (|A4p ]. In the /(i?) model, 
the background expansion is consistently determined by the action ([T]) once the form of f{R) is specified. In the 
present paper, without specifying the explicit form of f{R), we only adopted Eq. ^ in an asymptotic region. And 
we assumed the ACDM model as the background expansion of the Universe in the previous section. But we here 
consider possible uncertainties of the background expansion, by including the parameters wq and Wa- However, as 
will be shown in the below, the inclusion of the parameters Wq and Wa does not alter our result at the qualitative 
level. 

In the Fisher matrix analysis, we assume the galaxy sample of a survey with the number density Ng — 35 per 
arcmin.^, the mean redshift z,n = 0.9, and the total survey area, AA — 2 x 10'* square degrees. We also assumed 
the tomography with 4 redshift bins (see also Appendix) . Figure |4] is the result of the Fisher matrix analysis of 
the 9 parameters, n. A, wq, Wa, ^o-, ^b-, h, A, and Ug. Figure HI plots the 1-sigma contour in the n — \ plane, 
which is obtained by marginalizing the Fisher matrix over the other 7 parameters. The target values of n and A 
are shown in the panels. The other target parameters are wq = —1, Wa = 0, fio = 0.28, = 0.044, h = 0.7, 
Us = 0.95, and A which is set so that cts — 0.8. We take into account the Planck prior constraint of the expected 
errors Awo = 0.6, Awa = 1.9, Al^o = 0.01, Af^b = 0.0014, Ah = 0.01, Adg = 0.1 and An^ = 0.014 The left 
panel shows the result using the linear theory for the matter power spectrum of the range, 10 < ; < 10^. The right 
panel is the result with the nonlinear matter power spectrum of the range, 10 < Z < 3 x 10'^. In this figure, the solid 
curve corresponds to kc = 0.2 /iMpc~*jWhich was defined as the boundary between the general-relativity regime and 
the dispersion regime in the reference [67|. 
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FIG. 5: Left panel: The 1-sigma contour in the {kc — n) plane. The linear modeling for Pmass(fc, z) in the range of 10 < ? < 10^ 
is used. We assume the target modes kc = 0.1/iMpc~^ and n = 5, 4, 3, 2, 1, and 1/2, respectively, from the larger circle to 
smaller one. The other parameters are the same as those of Fig. [H Right panel: Same as the left panel, but with the nonlinear 
modeling for Pmass(fc, z) of the range of 10 < ? < 3 x 10^. 



Figure [5] is similar to Fig. |4l but the l-sigma contour in the kc — n plane. We assumed the target modes 
kc = 0.1/iMpc~^ and n = 5, 4, 3, 2, 1, and 1/2, from the larger circle to smaller one, respectively. The other 
parameters are the same as those of Fig. 21 Figure [5] shows the 1-sigma error on kc as a function of the target 
value of kc, where the other parameters are marginalized over. The left panels are the linear theory, while the right 
panels are the nonlinear model. The upper panels are the result of the Fisher matrix of the 9 parameters n, kc, 
wq, Wa, ^0, ^b, h. A, and rig. Then, the 1-sigma error Akc is evaluated by marginalizing the Fisher matrix over 
the 8 parameters n, wq, Wa, ^o, ^b, h, A, and Ug. The error of kc is the same order of kc for the cases n — 1/2 
and 1, but the error becomes larger as n becomes larger. The lower panels are the result of the Fisher matrix of the 
7 parameters, kc, n, flo, fif,, h. A, and Us, with fixing the background expansion to be that of the ACDM model. 
Thus, the inclusion of the parameters wq and Wa does not alter the result qualitatively. 

Figures HUB] show that the difference between the linear modeling and the nonlinear modeling is not very significant. 
We adopted the Peacock and Dodds formula [12] for the nonlinear modeling of the matter power spectrum, while 
the formula by Smith et al. 1931 has been used frequently [1^ |93 - [97| . However, the choice of the nonlinear formula 
doesn't alter our conclusion qualitatively. We have not taken the nonlinear effect from the Chameleon mechanism 
into account. The nonlinear modeling for the f{R) model has not been studied well for the general case of n. The 
effect of the nonlinear modeling might need further investigations. 



V. SUMMARY AND CONCLUSIONS 



In the present paper, we have investigated the linear growth rate of cosmological matter density perturbations in 
the viable f{R) model both numerically and analytically. We found that the growth rate in the scalar-tensor regime 
can be characterized by a simple analytic formula (jl5p . This is useful to understand the characteristic behavior of 
the growth index in the scalar-tensor regime. We also investigate a prospect of constraining the Compton wavelength 
scale of the f{R) model with a future weak lensing survey. This result shows that a constraint on kc of the same 
order of kc will be obtained for the model n = 1 and n = 1/2, though the constraint is weaker as n is larger. For 
kc ^ 1/iMpc^^, the constraint is very weak. This is because the weak lensing statistics is not very sensitive to the 
density perturbations on the smaller scales. 

Finally we mention about the effect of the late-time evolution of matter density perturbations in the /(i?) model on 
the spectral index. This effect causes the additional spectral index, which is evaluated by Aus = dlnDjja, k)/d\nk. 
The analytic formula of the additional spectral index is given by Starobinsky [6^ (see also [1^ IHI), Aug ~ 
(n/33- 5)/(6n4-4), which yields Aus = 0.11, 0.074, 0.047, 0.034, 0.027, and 0.022, for n = 1/2, 1, 2, 3, 4, 
and 5, respectively. Figure [7] plots our numerical result of An^ as a function of k assuming kc — O.l/iMpc"^. The 
numerical result approaches the analytic result at k ^ kc, but one can see the bump around the wavenumber kc, 
depending on n. Possibility of detecting of the spectral shape is interesting, but is out of the scope of this paper. 
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FIG. 6: The 1-sigma error on kc as a function of the target value of fee, where the other parameters are marginaUzed over. 
The left (right) panels use the linear (nonlinear) modeling for Pmass(fc, z) of the range of 10 < Z < 10^ (10 < / < 3 x 10^). 
In each panel, the curves assume the target parameter n = 5, 4, 3, 2, 1, and 1/2, from the top to the bottom, respectively. 
The other target parameters are the same as those of Fig. U The upper panels are the result of the 9 parameters, kc, n, 
wo, Wa, f2o, 0.b, h. A, and Ua. The lower panels are the result of the 7 parameters, kc, n, fio, ^b, h, A, and Us, where the 
background expansion is fixed as that of the ACDM model. 
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Appendix A: Modeling of weak lensing survey power spectrum 

We briefly review the Fisher matrix analysis for a weak lensing survey. The analysis in the present paper is almost 
the same as that of Ref. [l^, but the diff'erence is the modeling for the evolution of the matter density perturbations. 

As is described in Sec. IV, the signal of the weak lensing is determined by $ — ^' along the path of a light ray. 
Assuming the weak lensing tomography method [985 , the cosmic shear power spectrum for the i-th and j-th redshift 
bins is 

%)(0 = J^^ dxWMx))W,{z{x)) (^^) JP*-* (^k ^ ^,z(x)) , (Al) 
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-2-10 1 

logk[h Mpc"'] 



FIG. 7: The additional spectral index An^ as a function of wavenumber fc(/iMpc ^). Here we adopted kc ~ 0.1/iMpc ^ . The 
curves assume n — 1/2, 1, 2, 3, 4, and 5, from the top to the bottom, respectively. 



where P$_^(fc, z) is the power spectrum of $ — ^t, x is the comoving distance, Wi{z) is the weight factor of the i-th 
redshift bin, 

N^ Jn,^K(z,.z) dz' \ X(z')J 

where dN/dz{z) denotes the differential number count of galaxies with respect to redshift per unit solid angle, and 
Ni — J^^^ dz'{dN{z')/dz') is the total number of galaxies in the z-th redshift bin. From Eq. (pS)) . Eq. (jAip is written 
as 

P(,,)(0 = dxWMx)WMx)) (^^) ^mass [k ^ , (A3) 

where Praass{k,z) is the matter power spectrum, for which we adopted the Peacock and Dodds formula [o^l for the 
nonlinear modeling. This expression (jA3p is familiar as the weak shear power spectrum, but the modification of the 
gravity is involved in the evolution of the matter power spectrum Pmass(fc, z). 
In the present paper, we adopt the comoving distance 

~ Jo H{z') ^ Jo i7oVf^o(l + z'f + (1 - f^o)(l + z')3(i+"'"+"''')e-3'""^'/(Wy' 

which includes wq and Wa, the parameters of the equation of state of the dark energy w{z) = wq + Wa(l ~ a)- As 
mentioned in Sec. IV, the background expansion of the f{R) model is specified once the form of f{R) is given. The 
f{R) model, in the present paper, only assumes the form in an asymptotic region. Taking possible uncertainties of 
the background expansion, we include wq and Wa in the Fisher matrix analysis. However, this inclusion does not alter 
our result. 

The fisher matrix is estimated as 

I (ij){mn) 

where 9°' is a parameter of the theoretical modeling, and the covariance matrix is 

Cov(,,-)(_) (0 = [Pit) (0^(t) (0 + Pit) (O] > (A6) 
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with P°^^{1) = P(^ij){l) -\- Sijal / N i, where /sky is the fraction of the survey area, and CTe is the rms value of the intrinsic 
random eUipticity, which we take 0.22. In the Fisher matrix analysis, we assume the sample of galaxies of imaging 
survey modeled as 



dz 



Ng(3 



z exp 



(A7) 



with a — 0.5, /3 = 3, A'g = 35 per arcmin^, and zq is given by the relation, zo = z,„r((Q; + 1) / (3) /^{{a + 2) so that 
the mean redshift is z^ = 0.9. Wc assume the survey area, Av4 = 2 x 10^ square degrees, and the tomography with 
4 redshift bins. 
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